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ABSTRACT 

In this work we present constraints on different shapes of primordial non-Gaussianity 
using the Wilkinson Microwave Anisotropy Probe (WMAP) 7-year data and the spher- 
ical Mexican hat wavelet /„/ estimator including the linear term correction. In par- 
ticular we focus on the local, equilateral and orthogonal shapes. We first analyse the 
main statistical properties of the wavelet estimator and show the conditions to reach 
optimality. We include the linear term correction in our estimators and compare the 
estimates with the values already published using only the cubic term. The estimators 
are tested with realistic WMAP simulations with anisotropic noise and the WMAP 
KQ75 sky cut. The inclusion of the linear term correction shows a negligible improve- 
ment 1 per cent) in the error bar for any of the shapes considered. The results of 
this analysis show that, in the particular case of the wavelet estimator, the optimality 
for WMAP anisotropy levels is basically achieved with the mean subtraction and in 
practical terms there is no need of including a linear term once the mean has been 



f(/oc) 



39.0±21.4, f% q) 



-62.8± 154.0 and 



subtracted. Our best estimates are now: /J^ 

f\ — —159.8 ± 115.1. We have also computed the expected linear term correction 
for simulated Planck maps with anisotropic noise at 143 GHz following the Planck 
Sky Model and including a mask. The improvement achieved in this case for the local 
/„; error bar is also negligible (0.4 per cent). 
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1 INTRODUCTION 

In the rec ent years the spherical Mexi can hat wavelet 
(SMHW) (|Martmez-Gonzalez et~ai1 120021 ) has been used 
to construct a new type of estimator for the primor- 
dial non-Gaussianity in the CMB characterised by t he 
non-l i near c oupling parameter /„; l|Curto et alj l2009al lbl. 
|2010| . l2011al lbh. One of the particularities of the wavelet 
estimator as it has been traditionally presented in the 
literatur e compared with d i rect b ispect rum-based esti- 
mato rs ( Komatsu fc Spergell 200ll ; iKomatsu et all |200^ 
20031: iBabich et all 12004 iBabichl 120051: ICreminelli et al 



2006, 200 71: lYadav fc WandeltJ 120 08: Sm ith et al.l 2009: 
Eisner fc Wanderdl2009|; iLiguori et"a!Tl2010l; ISenatore et all 



20101: 



iSmidt et aJj 



Fcrgusson et al. 



2010al lbl: 



Komatsu et al.l l201ll : iFergusson fc Shellardl l201ll ) is the 
absence of a linear term. In the bispectrum-based estima- 
tors, the linear term plays a key role to achieve optimality 
in the cases where the rotational invariance of the CMB 



is broken because of different instrumental complexities 
such as anisotropic noise or pa rt ial sk y coverage (see for 
example ICreminelli et~ai1 |2006l. 120071 : lYadav fc Wandeltl 
l2010l : iFergusson fc Shellarcfeoill ). 

The computational difficulties related to the inversion of 
the covariance matrix present in the bispectrum estimator, 
especially in future data sets with higher l m ax as for exam- 
ple Planclfl, together with the unknown effect that different 
systematics from the instrument and background residuals 
might have on the estimates, motivated the search for new 
estimators based on different tools such as th e SMHW de- 
scrib ed in this paper, the binned bispectrum l|Bucher et al.l 
|2010| ). the general modal expan s ion and polyspectra 



estimation (IFergusson et all l2010bl; iFerg usson fc Shellardl 
2011), the needlets |Marinucci et alj 20081: Pietrobon et al 



20091: iRudiord et al 



2009 



Donzelli et alj I2012T ), the 



HEAL Pix wavelet Jcasaponsa et alj l2011al ). neural net- 



works (lCasaponsa et al.1 l2011bF or a Bay esian approach 
~alj|20ld : lElsner fc WandeltJl201ol) among others. 



( Eisner et 
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In a previous paper (|Curto et alj|2011a[ ). we described 
the main features of the wavelet estimator based on the cu- 
bic statistics constructed from the SMHW coefficient maps. 
Those cubic terms were written as a function of the non- 
linear coupling parameter /„; and the bispectrum of the 
primordial non-Gaussianity. In that paper we also showed 
that the power of the method to detect /„;, that is the 
variance of this parameter a 2 (f n i), matches that of the di- 
rect bispectrum-based estimators for ideal conditions (full 
sky and isotropic noise) and realistic conditions (partial 
sky coverage and anisotropic noise). The wavelet estima- 
tor variance was obtained in two different ways: through 
the Fisher matrix and by means of Monte Carlo (MC) sim- 
ulations, providing very similar results. A remarkable re- 
sult of these works is the fact that the wavelet estimator 
is, in practice, able to reach optimality on the /„; estima- 
tion without including any linear ter m correction. However , 
from several works (see for example ICreminelli et af1l2006l : 
iFergusson fc Shellardllioill ) it has been shown that in order 
to reach minimum variance, all the cubic estimators need 
a linear term corre ction. A recent work has solved this ap- 
parent controversy |Donzelli et al.ll2012l ) by showing that in 
anisotropy conditions, the linear term correction 
is nearly equivalent to the mean subtraction performed for 
each wavelet coefficient map in the wavelet estimator. 

In this paper we re-examine the main statistical prop- 
erties of the wavelet estimator and show the conditions to 
reach optimality. We compute the linear term correction for 
the local, equilateral and orthogonal /„; shapes. In particu- 
lar we see that the linear term correction for the local case 
provides a 1 per cent reduction in the error bars (in agree- 
ment with iDonzelli et~aT1l2012l , for the SMHW) while the 
correction for the other shapes is even smaller. Section [2] in- 
troduces the SMHW estimator, its variance and its linear 
correction. In Section [3] the estimator with its linear correc- 
tion is applied to WMAP 7-year data for the local, equilat- 
eral and orthogonal shape. In Section[4]we explore the linear 
correction on Planck simulations at 143 GHz for the local 
shape, and in Section [5] the conclusions are presented. 



2 THE WAVELET APPROACH 

In this Section we present an approach for the f n i estimator 
based on the statistical properties of the cubic terms of the 
SMHW coefficients averaged over the sky. In this case we 
exploit the property of the SMHW wavelet that performs a 
strong decorrelation of the data at distances larger than the 
wavelet resolution. The expected value of the cubic terms 
in the sky are obtained from the sum of a large number of 
almost independent elements and therefore its distribution 
will be close to Gaussian by the central limit theorem. We 
will first review the SMHW and its decorrelation properties 
and then we will construct the wavelet estimator based on 
those properties including the linear term correction. 




.0 10.0 
B (deg) 

Figure 1. The normalised correlation functions for different 
wavelet coefficient maps are plotted. From left to right, the curves 
correspond to maps convolved with a SMHW of scale R = 0.27, 
0.64, 1.50, 3.54 and 8.33 degrees, respectively. 



2.1 The SMHW coefficients and their correlation 

Detailed information about the spherical Mexican Hat 
wavelet (SMHW) and a (non-complete) list of applica- 



lAntoine & Vandergh 


evnstl ( 


1998); Martinez-Gonzalez et al. 


(2002); Cavon et al.1 (2003 


); IVielva et al. 


d200eh: 


Viclva 


(2007); McEwen et al. (2007); Martinez-Gonzalez! ( 


2008); 



IZhang et al.l (|201lD ; lYu et al.l (|2012D . 

Given a function /(n) defined at a position n on 
the sphere and a continuous wavelet family on that space 
*(n; b, R), we define the continuous wavelet transform as 



w(R;b)= / dn/(n)*(n;b,7?) 



(1) 



where b is the position on the sky at which the wavelet 
coefficient is evaluated, R is the scale of the wavelet and 
* s (6»; R) = *(n(e, <f>); 0, R) is given by 



9 S {9;R) 
where 

N(R) = R ^1 - 
and 

y = 2 tan (~^J 



2tyN(R) 



1 + 



2 - 



ay 



-y 2 /{2R 2 ) 



(2) 



R 2 R 4 



1/2 



(3) 



(4) 



Considering a set of different angular scales {Ri} we de- 
fine a third order st atistic depending on three scales {i, j, k} 
l|Curto et al.ll2009bl ) 



Qijk 



1 



47T OiO ,'CTfc 



dnw(Ri, n)w(Rj, n)w(R k , n) 



(5) 



where Oi is the dispersion of the wavelet coefficient map 
w(Ri,n). In the particular case of Ro = 0, w(Ro, n) = / (n). 
For a particular pixelization on the sphere, Eq. [5] can be 
written as: 



http:/ /map. gsfc.nasa.gov/ 
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qijk = 



N 



N pix -1 

y> Wj(p)wj(p)w k (p) 

p=0 



(6) 



where N p i x is the total number of pixels of the map, Nijk is 
the number of pixels available after combining the extended 
masks corresponding to the three scales Ri, Rj and Rk and 
u>i(p) = w(Ri,p) — (w(Ri)} is the wavelet coefficient in the 
pixel p evaluated at the scale Ri after subtracting the mean 
value over the wavelet coefficient map outside its extended 
mask. 

Using the properties of the wavelet, we may write the 
wavelet transform o f the temperature map in the next form 
(|Curto et al.ll2011al ) 



w(Ri, n) = y air, 



.Ui(Ri)Y em (n) 



(7) 



Using the isotropic properties of the CMB and the properties 
of the wavelet, we can obtain the angular coefficient correla- 
tion dj (9) between any pair of pixels n and n' separated by 
an angular distance nn' = cos(#) and for two angular scales 
Ri and Rj 



= {w(Ri, n)w(Rj,n)) = 
C t ut{Ri)u>t(Rj)Pt (cos(0)) 



(8) 



where uje(R) is the window function of the wavelet at a 
scale R and it is given by the h armonic transform of the 
mother wavelet of th e SMHW |Martmez-Gonzalez et all 
120021 ; ISanz et alj|2006f ). The dispersion of the wavelet co- 
efficients at scale Ri is simply given by Oi = Cu{9 = O) 1 ^ 2 . 

In Fig. [T] we show the correlation of the wavelet coef- 
ficients as a function of the angular distance 9 for several 
values of the resolution scale R. As can be seen, the SMHW 
produce an effective decorrelation of the signal at angular 
distances above the resolution scale R. 



2.2 The wavelet estimator 

Considering the strong decorrelation produced by the con- 
volution of the SMHW on the temperature anisotropics, we 
can now apply the central limit theorem to the cubic statis- 
tics defined in Eq. Since the average value is calculated 
from the sum of a very large number of almost independent 
elements (of the order of the number of pixels in the sphere 
with size that of the resolution scale R), then its distribu- 
tion should be very close to a Gaussian. This is actually 
seen in Fig. O where the distribution of the cubic terms for 
different SMHW scales are shown. These distributions have 
been obtained from MC simulations of Gaussian tempera- 
ture anisotropies. 

The previous results indicate that, for Gaussian tem- 
perature anisotropies, a good representation of the n-point 
distribution of the quantities qijk can be given in terms of a 
multinormal distribution. Allowing now for the presence of 
weak non-Gaussianity for the temperature anisotropies (e.g. 
an amplitude for the primordial non-Gaussianity consistent 
with WMAP data) one can use the next likelihood for the 
f„i parameter: 



L(f nl ) oc e-* 2 <'»'>' 2 , 
where x 2 (fm) is given by 



(9) 



X 2 (fnl 



= 

ijk,rst 



bs 
ijk 



fnlOlijk)C iikrst {q rat 



fnlCtrst), (10) 



where are the cubic statistics corresponding to the ob- 
served data, ctijk = {Qijk) f„i=i an d Cijk,rst is the covariance 
matrix of the cubic statistics. A further test to check that 
the qijk are normally distributed can be done by consider- 
ing the property that Ax 2 (/ni) = X 2 (.fm) - XmiMm) is a 
X 2 distribution with one degree of freedom. In particular, 
Ax 2 (/n0 = 1(4) should provide the 1(2) -sigma or 68 per 
cent(95 per cent) confidence intervals for the /„; parameter. 
Using MC simulations we have checked that this is the case 
for the qijk statistics. 

After straightforward calculation, it can be easily seen 
that the /„; estimator in this case is given by 



fnl — 



^2 a ijkC i . k rst q rs . 

X] a 'i k ^-'ijk,rst ars 



(11) 



while that the variance of the /„; parameter in Eq. (|11[) is 
given by 



J" anogL{f nl ) \ 

I 8 '"' //„,=/„ 



(12) 



l f W \ Eijk, rst ^C-f rat a rst ■ 

L S/ " ! J fnl=fnl 

This estimator h as already been show n to be nearly optimal 
on WMAP data (|Curto et al.ll2009d lbl l20ldl2011al lbl) with- 
out the need of subtracting a ny linear term. However, as 
stated in iDonzelli etall (|2012t) . from all the possible cubic 
combinations of three Gaussian variables, the Wick poly- 
nomials are shown to have minimum variance. This implies 
that in order to have a strictly speaking minimum variance 
estimator, a linear term correction needs to be included. In 
fact the linear term subtraction is equivale nt to the mean 
subtr action at each wavelet coefficient map l|Donzelli et al.l 
l2012h for low levels of anisotro py. This is indeed th e proce- 
dure t hat has been followed in ICurto et ail (|2009al lbl |2010| . 
l2011al fbT) and it explains the competitive results obtained 
just by subtracting the mean using the estimator in Eq. 1111 
The linear term correction for the wavelet estimator can 
be written as: 



r(total) r(cubic) ?(Zinear) 

J nl •> nl J nl 



where /^j" 6 ^ is given by Eq. Illl and 

r(linear) 2 / f \ \ * r 1 ~^ C^) 

Jnl — a F\Jnl) / J &ijkLsijk,rst C lrst 

ijk,rst 



(13) 



(14) 



with 



lljk 



1 

Nijk 



E 



Wi{p) Wj(p) \ w k {p) 
Oi Oi I 



Ok 



Wi (p) W k (p) \Wj(p) I Wj (p) Wk (p) \ Wi (p) 



Ok 



+ 



Ok 



(15) 



In the next Sections we apply the wavelet estimator to 
WMAP 7-year data as well as to Planck simulations and 
compare the results obtained with and without the linear 
term correction. 
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3 APPLICATION TO WMAP V+W DATA 

We have computed the linear term correction to the cubic 
wavelet f n i estimator for the three shapes with a relevant 
interest in many inflationary models: t he local, equilateral 
and orthogonal shap e s (see for example |Bartolo et al, 1 12004 
ISenatore et al.l |2010| ; iKomatsu et al.l 12011 ) . The estimator 
can be easily applied to other bispectra with a separable 
shape (|Curto et alj|2011ah . Results taking into account only 
the cu bic contribution are already published IjCurto et al.l 
l2011bh . 

We have selected the same set of 15 angula r scales from 
Ro = arcmin to R14 — 767.3 arcmin used in lCurto et al.l 
l|2011bf ). We have considered the V+W WMAP data op- 
timally weighted by the Nhits maps per radiometer in or- 
der to maximise the signal-to-noise ratio. We also consider 
the same WMAP KQ75 mask and its extended masks for 
each wavelet angular scale. The cubic covariance matrix 
has been computed using 10000 Gaussian simulations. A 
principal component analysis has been performed in order 
to avoid contamination from the lowest noisy eigenvalues 
of t his covariance matr ix without losing non- Gaussian sig- 
nal (|Curto et al.l l2011al ). The two point correlation matri- 
ces needed for the linear term correction have also been es- 
timated with 64000 Gaussian simulations. This number of 
simulations is needed in order to achieve the required preci- 
sion in the estimation of the correlation matrices. 



We have applied the estimator to one set of 10000 Gaus- 
sian maps and the WMAP data. The results are presented 
in Fig. [3] for the three considered shapes. In the left pan- 
els, the red histograms correspond to the best-fitting f n i 
values obtained with the cubic estimator and the black his- 
tograms correspond to the best-fitting /„; values after the 
linear term correction. The vertical lines correspond to the 
actual WMAP data values estimated with the cubic esti- 
mator (red) and the linearly corrected estimator (black). 
In the right panels, we compare the best-fitting f„i values 
for the same set of Gaussian simulations. Note that both 
/iz" 6lC ^ ano - fni"* are highly correlated and the deviations 
are not significant. Finally in Tables [T] [2] and [3] the previous 
results are summarised. We present the WMAP 7-year /„/ 
best-fitting values for the cubic estimator f^ l ublc ' > , the linear 
estimator /^; lnear ', and f^° tal ^ for the clean and raw (un- 
clean ed) maps. The Fishe r f n i error bar as described in Eq. 
25 of ICurto et al.l (|2011al ) is also provided. For each case, 
we observe a small reduction of the error bars when the lin- 
ear term is included. The largest correction is introduced in 
the local shape, where a is reduced from tj(fni) ~ 21.6 to 
c(.fn i) = 21.4 (i.e. a red uction of 1 per cent, in agreement 
with iDonzelli et ai1l2012l ). The correction for the other two 
cases, equilateral and orthogonal is also negligible (about 
0.2 per ce nt and 0.1 per cent resp ectively). This is in agree- 
ment with lCreminelli et al.l (|2006l ) for the equilateral shape, 
where the standard deviations of f n i without the linear term 
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-400 -200 200 400 -400 -200 200 400 

WMAP best — fitting equilateral f n] equilateral cubic f n . 




WMAP best — fitting orthogonal f n orthogonal cubic f n . 



Figure 3. WMAP 7-year data best-fitting /„; values using the cubic estimator (Eq. [TT} and the estimator with the linear correction (Eq. 
1 1 3 D for the local (top), equilateral (middle) and orthogonal (bottom) shapes. In the left panels the histogram of the best-fitting /„; values 
with (dashed dark line) and without (solid red line) the linear term correction for each simulation are plotted. Vertical lines correspond 
to the values obtained with WMAP data. The right pannels show the corresponding correlation between the same estimates. 



were found closer to the lower Fisher limit than in the local 
shape, suggesting a less important contribution of the linear 
term correction. 

Our best-fitting values, computing o(f n i) with 10000 
Gaussian simulations to characterize the errors, are pre- 
sented below for the three shapes. 

Local form results^: 

• fni Ubic) =38.9 ±21.6 

3 The average (<}ijk)f n i is obtained using 1000 non-Gaussian 
simulations of the local shape gene rated by the procedure de- 
scribed in jElsner fc Wandehj 2009} and publicly available at 
http: / / planck. mpa-garching. mpg.de / cmb / fnl-simulations / 

The best estimates of the local shape presented in lCurto et al.l 
ll2011bh are f^ ubic) = 32.5 ± 22.5. Note that in that work, a 
perturbative approach is considered to simulate the non-Gaussian 
simulations used to compute (qijk)f l - The different approaches 
to simulate the non-Gaussianity and the statistical errors due 
to the finite number of non-Gaussian simulations explains the 
small differences between the error bars presented here and in 
that reference. 



• flT al) = 39 -° ± 21 - 4 

Equilateral form results: 

• fn"^ = ~ 53 - 3 ± 154 - 3 

• f%° iaV) = -62.8 ± 154.0 
Orthogonal form results: 

• /£T WC) = -155-1 ± 115-1 

• fn° tal) = -159.8 ±115.1 

In order to check that the estimator has already reached op- 
timality with the considered scales for the three shapes, we 
have computed a(fni) for different subsets of scales (Fig. [4}. 
We compare the /„; error bars for different minimum angular 
scales Rmin- To find the equivalent multipole I range cor - 
responding to each R m in see Fig. 5 of ICurto et al.l (|2011al ). 
The three shapes reach minimum variance for Rmin ~ 
arcmin. 

The error bar of the equilateral and orthogonal shapes 
are also similar to the values obtained with the direct bispec- 
trum estimator where a(f n i) = 140 for the e quilateral shape 
and a (fni) = 104 for the orthogonal shape (jKomatsu et al.l 
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Table 1. Constraints on the / n ; parameter for the local shape with and without the linear term correction. From left to right, the 
best-fitting values for the clean and the raw data maps, the mean, dispersion, 16, 84, 2.5 and 97.5 per cent quantiles respectively of the 
f n l distribution obtained with 10000 Gaussian maps. The Fisher error bar obtained for this shape is <Tp(f n i) = 21.6. 



case 


/.(clean data) 


{.(raw data) 
Jul 


(fnl) 


v(fnl) 


X\6 


^84 


X2.5 


^97.5 


cubic 


38.9 


20.8 


0.6 


21.6 


-21.1 


21.9 


-42.7 


41.8 


linear 


-0.1 


-0.0 


0.0 


3.1 


-3.1 


3.2 


-6.1 


6.2 


cubic - linear 


39.0 


20.8 


0.7 


21.4 


-21.0 


22.7 


-42.6 


41.3 



Table 2. Constraints on the /„; parameter for the equilateral shape with and without the linear term correction. From left to right, the 
best-fitting values for the clean and the raw data maps, the mean, dispersion, 16, 84, 2.5 and 97.5 per cent quantiles respectively of the 
/„; distribution obtained with 10000 Gaussian maps. The Fisher error bar obtained for this shape is crp(f nl ) = 144.5. 



case 


(.(clean data) 
Jnl 


p(raw data) 
Jnl 


(fnl) 


"(fnl) 




^"84 


X2.5 


^97.5 


cubic 


-53.3 


28.1 


-1.6 


154.3 


-155.9 


151.5 


-302.4 


302.2 


linear 


9.5 


13.7 


-0.3 


23.0 


-23.6 


22.6 


-47.8 


46.7 


cubic - linear 


-62.8 


14.4 


-1.3 


154.0 


-156.4 


150.3 


-304.5 


300.3 



Table 3. Constraints on the /„; parameter for the orthogonal shape with and without the linear term correction. From left to right, the 
best-fitting values for the clean and the raw data maps, the mean, dispersion, 16, 84, 2.5 and 97.5 per cent quantiles respectively of the 
f n l distribution obtained with 10000 Gaussian maps. The Fisher error bar obtained for this shape is crp(f n i) = 106.3. 



case 


Aclean data) 
Jnl 


f.(raw data) 
Jnl 


(fnl) 


"(fnl) 


X\6 


^84 


X2.5 


^97.5 


cubic 


-155.1 


-119.4 


0.2 


115.1 


-113.4 


115.2 


-230.4 


225.6 


linear 


4.7 


4.8 


0.2 


18.5 


-17.8 


18.0 


-36.5 


36.7 


cubic - linear 


-159.8 


-124.2 


0.0 


115.1 


-113.2 


115.0 


-228.4 


222.8 





cubic 




+ cubic - 


linear 





125 - 


b 








120 - 


cn 







115- 



5 10 



15 20 



Figure 4. The /„; error bars computed using different K m in angular scales for the local (left), equilateral (middle) and orthogonal 
(right) shapes. The diamonds corresponds to f™ blc and the crosses to f t ^, tal ■ 



2011). The slightly larger values (~ 9 per cent) obtained 
from the dispersion of the /„( distribution corresponding to 
10000 Gaussian simulations, a(f n i) = 154 and u(f n i) — 115 
respectively, are likely due to differences in the perturbative 
approach us ed to simulate the non-Gaussian signal of these 
two shapes |Curto et al.ll2011bh or the statistical errors due 
to the finite number of non-Gaussian simulations. 



4 APPLICATION TO PLANCK SIMULATIONS 

We have computed the linear term correction to the cubic 
wavelet f n i estimator for the local shape using Planck simu- 
lations in order to forecast the amplitude of this correction 
on future Planck analyses. We do not consider the two other 
shapes (equilateral and orthogonal). From the results of pre- 



vious sections, we expect the correction for those cases to be 
even smaller. 

For this analysis we have considered a new set of angular 
scales that better suits the range of angular multipoles which 
are cosmic variance dominated (^ ma x ~ 2000). The list of 
angular scales is Rq = 0, Ri = 1.3, R2 = 2.1, R3 = 3.4, 
R 4 = 5.4, R 5 = 8.7, R 6 = 13.9, R 7 = 22.3, R s = 35.6, R 9 = 
57.0, Rio = 91.2, Ru = 146.0, R12 = 233.5, R13 = 373.6, 
i?i4 = 597.7 and R15 = 956.3 arcmin. As a representative 
mask, we have used the available WMAP KQ75 mask (75 
per cent of the sky). We have simulated the Planck 143 GHz 
channel using a fiducial CMB power spectrum that best fits 
WMAP 7-year data, l m ax = 2048 and a Gaussian beam 
with FWHM — 7.1 arcmin. The noise has been generated 
using an anisotropic Nhits map computed from the scanning 
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Table 4. Expected constraints on the /„; parameter for the local 
shape with and without the linear term correction for the 143 GHz 
Planck channel. From left to right, the mean, dispersion, 16, 84, 
2.5 and 97.5 per cent quantilcs respectively of the / ni distribution 
obtained with 1000 Gaussian maps. The Fisher error bar obtained 
for this case is a"i?(/ n ;) = 7.5. 



case 


(fnl) 




^16 


^84 


^2.5 


-^97.5 


cubic 


0.4 


7.98 


-7.2 


8.6 


-15.0 


17.2 


linear 


0.0 


1.43 


-1.4 


1.5 


-2.7 


2.8 


cubic - linear 


0.4 


7.95 


-7.4 


8.2 


-15.4 


16.7 



strategy of the Planck Sky Modefl ijDelabrouille et alJl20l3 ) 
and the noise sensitivity per pixel provided in the Planck 
Bluebooljf] (using an average noise sensitivity for 14 months 
of (Tnoise = 2.2 /xK/K in a square pixel whose size is the 
FWHM extent of the beam). 

The cubic covariance matrix and the linear correlation 
matrices needed for the f n i estimator in Eq. (|13[) have been 
computed using two independent sets of 10000 Planck Gaus- 
sian simulations. The results corresponding to the analy- 
sis of an additional set of 1000 Gaussian maps are pre- 
sented in Fig. [S] Note that for this simulated Planck level 
of anisotropy, /„; Ui>!C ' and /„j 0< are also highly correlated. 
Finally in Table U the properties of the previous histograms 
are summarised. In particular, we see that using the cubic 
estimator a(f n i) = 7.98 and the linear term contribution 
reduces this error bar to a(f n i) = 7.95 (i.e. a negligible cor- 
rection lower than 0.4 per cent). 



5 CONCLUSIONS 

In this paper we have performed a comprehensive study of 
the f n i wavelet estimator. We have considered the main sta- 
tistical assumptions to derive the estimator in terms of the 
cubic quantities qi and showed the conditions to reach op- 
timality. We have found that the SMHW wavelet produces 
an important effective decorrelation of the signal at angular 
distances above the angular resolution R. This means that 
the cubic quantities are nearly Gaussian owing to the cen- 
tral limit theorem and using this property we have found 
an expression for the likelihood of the f„i parameter in 
terms of the cubic statistics. We have also included a linear 
ter m correction fo l lowin g the Wick polynomials introduced 
by iDonzelli et al.1 l|2012ft . In particular, we have confirmed 
that the linear term correction is basically achieved through 
the mean subtraction that we carry out on the wavelet co- 
efficient maps for each angular scale (|Curto et alj|2011al lbh. 
We find that, in this case, the linear term correction only re- 
duces the error bars about 1 per cent for the local case using 
WMAP data. This correction is even smaller for the equilat- 
eral and orthogonal cases (0.2 per cent and 0.1 per cent re- 
spectively). The results presented in this paper are in agree- 
ment with the optimal results obtained with the wavelet es- 
timator alr eady published whe re the mean su btraction was 
performed (|Curto et alj|2009al lbl l20ld l2011al lbh. Therefore, 

4 http:/ /www. ape. univ-paris7.fr/~delabrou/PSM/psm. html 

5 The Planck Bluebook is available for download in the web: 
http: / /www.rssd.esa.int/index.php?project=Planck 



we conclude that the contribution of the linear term is neg- 
ligible (^ 1 per cent) for the SMHW estimator for the three 
considered shapes. We have also explored the linear term 
correction for Planck simulations at the 143 GHz channel. 
Our results indicate that the correction for the local shape 
is lower than 0.4 per cent considering the expected levels 
of noise anisotropy for this channel and the WMAP KQ75 
mask. From the results on WMAP data, we expect the cor- 
rection for the equilateral and orthogonal shapes to be even 
smaller. 
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